Heat capacity of liquids: an approach from the sohd phase 
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We calculate the energy and heat capacity of a liquid on the basis of its elastic properties and 
vibrational states. The experimental decrease of liquid heat capacity with temperature is attributed 
to the increasing loss of two transverse modes with frequency uo < 1/r, where r is liquid relaxation 
time. In a simple model, liquid heat capacity is related to viscosity and is compared with the 
experimental data of mercury. We also calculate the vibrational energy of a quantum liquid, and 
show that transverse phonons can not be excited in the low-temperature limit. Finally, we discuss 
the implications of the proposed approach to liquids for the problem of glass transition. 



INTRODUCTION 

Heat capacity is one of the most important physi- 
cal properties of a system, because it holds information 
about its degrees of freedom. Heat capacity is well under- 
stood in gases and solids. Liquids, on the other hand, re- 
main an exception, with the result that their heat capac- 
ity is often barely discussed in classic statistical physics 
text books fl| or texts dedicated to liquids [2,l3|,|j,l5|,yjLZ|- 
The general difficulty in calculating the energy and other 
thermodynamic quantities of a liquid is that, unlike a 
gas, it has strong interactions, but at the same time does 
not have small atomic vibrations as a solid [1.]. Because 
the inter-molecular interactions in a liquid are strong, the 
properties of this interaction, and, therefore, properties 
of the liquid, are strongly system-dependent pj]. Conse- 
quently, it is argued that no general recipe for calculating 
the energy of a liquid can exist [l| . 

Experimental data from many elemental liquids show 
that their constant- volume heat capacity, C„, decreases 
with temperature, from about 3A'' at the melting point 
to (2 - 2.5)A^ at high temperature [i,[i|. Here, N is the 
number of particles and ks — I- The decrease of Cy is 
also seen in many complex liquids [lO|. This behaviour 
is not understood in a consistent framework similar to 
those that exist for solid and gas phases. 

Most calculations of liquid energy, E, approach a liq- 
uid from the gas phase, and calculate the potential en- 
ergy in addition to the gas kinetic energy, giving E — 
K + U, where K is the kinetic energy of an ideal gas 
{K = 3NT/2) and U is the potential energy of inter- 
atomic interactions [l|, 12, S l5| . The simplest treatment 
assumes the case of a dilute system with only pair interac- 
tions. Assuming that the interactions are weak (except at 
short separations), the high-temperature expansion gives 
the system described by the Van Der Waals equation []| . 
While this approach can describe a dense gas or a liq- 
uid close to the critical point, it is not adequate for real 
liquids. For example, C„ of the Van Der Waals is equal 
to that of an ideal gas [l|], in contrast with the experi- 
mental results. The improved calculations of U employ 
higher-order (three- and four-) particle correlations [5(. 
The results, however, are not straightforward to use for 



numerical estimates i5| and require the knowledge of cor- 
relation functions as well as the properties of interatomic 
interactions. 

A distinctly different, and less common, approach 
to discuss the energy and heat capacity of a liquid 
is to include the strong interactions from the outset, 
by approaching a liquid from the solid phase. The 
temperature-dependent term of the energy of a solid is 
given by the phonon energy. In an isotropic solid (glass) 
all vibrations can be represented by one longitudinal and 
two transverse waves. In classical case, this gives the heat 
capacity of C„ = 3A^, the Dulong-Petit result. The dis- 
cussion then proceeds to establish how this result needs 
to be modified on transition from the solid to the liquid 
phase. Therefore, this approach accounts for the strong 
interactions in a liquid from the outset. In this sense, it 
is opposite to the approach that starts with the energy of 
an ideal gas and introduces interactions as a correction. 

In this approach, Brillouin modified the solid-like re- 
sult, C'v — 3iV, by assuming that a liquid does not 
support transverse waves, and encountered an interest- 
ing contradiction jll|. The total energy of a solid is 
E = {NT/2 + NT/2) + 2{NT/2 -f NT/2) = 3NT, where 
the first and the second term give the energy of longitu- 
dinal and transverse vibrations, respectively, and NT/2 
is the mean potential or kinetic energy. If a liquid loses 
two transverse vibrations, their potential terms vanish, 
but kinetic terms remain, and the liquid energy becomes 
E = {NT/2 + NT/2) + 2{NT/2) = 2NT , giving C„ = 2N 
,11]. This contradicted the experimental result that liq- 
uid heat capacity at the melting point is about the same 
as that of crystals, Cy — 3N. One way to resolve the con- 
tradiction is to assume that a liquid consists of small crys- 
talline domains that support transverse waves in some 
directions, giving C*,, = 3iV, but not in others, enabling 
a liquid to flow illl|. 

Cy was also studied in molecular dynamics simulations, 
with liquid Al as a case study |12l |. Similar to the exper- 
iment, the decrease of heat capacity with temperature 
was observed, and was interpreted as the progressive loss 
of liquid shear resistance [l^j . 

The idea that a liquid loses two transverse vibrational 
modes was also used by Landau to calculate the energy 



of phonon excitations in a quantum liquid at very low 
temperature. This gives the result that the temperature- 
dependent energy term of a liquid is three times smaller 
than that in a harmonic solid, reflecting the fact that 
only one longitudinal mode is preserved [l| . 

More recently, liquid heat capacity has been discussed 
on the basis of alternative mechanisms, including the con- 
sideration of liquid potential energy landscape with in- 
tervalley motions [9|]. 

In this paper, we propose that energy and heat ca- 
pacity of liquids can be understood on the basis of their 
solid- like elastic properties. In particular, we discuss how 
the idea of relaxation time r can be used to describe 
liquid vibrational states by approaching them from the 
solid phase. The experimentally observed decrease of liq- 
uid heat capacity with temperature is attributed to the 
increasing loss of two transverse waves with frequency 
u) < 1/r. In a simple model, we relate liquid heat ca- 
pacity and viscosity, and compare this relation with the 
experimental data of mercury. In addition, we calculate 
the vibrational energy of a quantum liquid, and show 
that transverse phonons can not be excited in the low- 
temperature limit. Finally, we discuss the implications 
of the proposed approach to liquids for the problem of 
glass transition. 



EXPERIMENTAL EVIDENCE FOR 
PROPAGATING MODES IN LIQUIDS 

The existence of dispersion relations and the ability to 
support propagating shear modes has been traditionally 
assigned to solids. Relatively recently, experiments have 
shown that these features also exist in liquids, with the 
evidence coming primarily from inelastic X-ray, neutron 
and Brillouin scattering experiments. In this section, we 
briefly recall some of these results. 

Early experiments detected the presence of propagat- 
ing modes in liquids, in the form of "dispersion curves" at 
high frequency, at the temperature around melting [l3|. 
Later, the measurements were extended to temperatures 
considerably above the melting point, and confirmed the 
existence of collective excitations in many liquids (see, 



for example, Refs. [IJ, |15|, |16|, ll7|). As reviewed re- 



cently, it is now well established that dynamics in liq- 
uids shows solid-like character, in that liquids can sustain 
high-frequency propagating modes down to wavelengths 
on the atomic scale, with solid-like dispersion relations 

The solid-like ability of a liquid to support high- 
frequency shear waves considerably above the melting 
point was directly observed some time ago [l7[. In ad- 
dition, the ability to support shear waves has been in- 
ferred from the widely observed "positive dispersion" , 
the increase of high-frequency sound velocity over that 
expected on the basis of hydrodynamics llj, |l5|, llq . 



The general explanation of this effect is that at high 
frequency, a liquid responds elastically to shear stress, 
developing an ability to support solid-like shear vibra- 
tions, in addition to longitudinal ones. As a result, the 
propagating speed starts approaching that in a solid, i.e. 
increases. Initially, these observations were made close to 
the melting point, and were attributed to the proximity 
of solid phase. However, later experiments found posi- 
tive dispersion at temperatures considerably above the 
melting point, confirming the general idea that liquids 
support solid-like shear waves at high frequencies that 
extend down to microscopic wavelengths [l6|. 



VIEW OF A LIQUID FROM THE SOLID PHASE: 
RELAXATION TIME 



Relaxation time of a liquid was phcnomcnologically in- 
troduced by Maxwell in the viscoelastic picture of fiow 
[iSj . This picture is based on the assumption that de- 
formation of a liquid in response to stress can be inter- 
polated as the sum of elastic and viscous terms, giv- 
^Pxy, where P^y, G and 77 are 



-g1& 



1 dP^y 



— o dy U df ■ 11 -^w - -^yj - / 

shear stress, shear modulus and viscosity, respectively, 
and the first and the second terms represent elastic and 
viscous response, respectively. When external motion 
stops, Vx — 0, stress relaxes as P = Pq exp(— f/r), where 
T = ij/G is the Maxwell shear relaxation time. 

Frenkel offered microscopic interpretation of r as the 
time between elementary rearrangement processes in a 
liquid [2|. He started the discussion with the contra- 
diction of essentially the same origin as encountered by 
Brillouin: on one hand, in order to explain the solid-like 
value of heat capacity in liquids at the melting point, 
Cv = 3N, thermal motion in liquids should be considered 
to be solid-like, i.e. vibration around fixed positions. On 
the other hand, this picture contradicts fluidity of liq- 
uids. Frenkel reconciled the contradiction by introducing 
T as the time between consecutive atomic jumps from 
one equilibrium position to another in a liquid. If t is 
large compared with the period of atomic vibrations, a 
liquid is characterized by vibrational states as a solid, 
giving Cy = 3N. If, at the same time, t is small com- 
pared with time t during which an external force acts 
on a liquid, usual liquid flow takes place. If, on the 
hand, r is large compared with t, a liquid responds to 
an external force elastically as a solid. Hence the elastic 
and viscous regime correspond to lut ^ I and wr -C 1, 
respectively, where uj is the frequency of external force. 
This can be further illustrated by considering an external 
force F = Acx-p{iujt), for which Maxwell's interpolation 
gives ^ = i(l -|- iujT)F, hence ujt ^ 1 and wr <C 1 
corresponds to elastic and viscous response, respectively. 

In this picture, Frenkel arrived at an important con- 
clusion about how the liquid state can be viewed and 
discussed. He proposed that a liquid should be viewed 



as an essentially elastic medium, in which elasticity is 
"masked" by fluidity at times larger than r (see Ref . [2^ ) . 
In this approach, Frenkel proposed to consider its vi- 
brational states using the following continuity argument. 
The solid vibrational energy is calculated for an isotropic 
solid, glass, and is expressed in terms of normal modes 
[l|. Frenkel noted that in his picture of liquid flow, no 
qualitative distinction can be made between a glass and a 
liquid: a solid glass is different from a liquid by the value 
of r only (see Ref. [^). Because the transition between 
a liquid and a glass is continuous, vibrational states of a 
liquid just above the glass transition temperature should 
not substantially change from those of a glass [4I . This 
argument points to the continuity of vibrational states 
between a solid and a liquid, and offers an important 
insight into a liquid from the solid phase. 



RELAXATION TIME AND VIBRATIONAL 
STATES IN A LIQUID 

At the phenomenological level, two ways to incorpo- 
rate T into the existing theories have been proposed. The 
first one starts from a solid phase, and modifies elasticity 
equations. This gives rise to frequency-dependent elastic 
moduli, which are interpolated between zero at wr <^ 1 
and their solid-state value at ujt 3> 1 [3 • Another way is 
to start from a liquid phase, and to modify the hydrody- 
namic equations, by either directly introducing Maxwell 
viscoelastic term in the Navier-Stokes equation Q] , or by 
introducing the exponential decay, with t as decay time, 
of the memory function of current correlations [J] . 

As reviewed recently, the microscopic understanding 
of liquid dynamics is still lacking [16|. Because elas- 
ticity and hydrodynamics describe a system at length 
scale that is large compared to interatomic separations, 
these approaches do not address the microscopic effects 
of relaxation: elementary relaxation processes in a liq- 
uid are local atomic jumps [l9|. Other disadvantages of 
the phenomenological description have been discussed in 
Ref. [l6|. An open question remains as to how can the 
liquid ability to sustain (lose) transverse modes with fre- 
quency uj > l/r (a; < 1/t) be explained at the atomistic 
level, similar to the normal-mode analysis in a solid? We 
propose that one way to understand the absence of propa- 
gating shear modes with frequency w < l/r in a liquid at 
the atomistic level is to note that local atomic rearrange- 
ment processes effectively break some of the interactions 
between neighbouring atoms and stop wave propagation. 

A local structural rearrangement involves a jump of 
an atom from its cage, with subsequent relaxation of the 
local structure [19|. We call these rearrangements local 
relaxation events (LREs). Lets consider how periodic 
transverse motion of atom i in a liquid affects the neigh- 
bouring atom j. Lets assume that atom i is vibrated 
with period tq in the direction perpendicular to the line 



connecting i and j, exciting a shear wave, and consider 
atom j is involved in a LRE with a period t. If atom i 
vibrates with period tq ^ r, its displacement does not 
affect atom j, because during time Tq atom j has inde- 
pendently jumped many times due to thermally-excited 
LRE. In this case, the wave does not propagate. If, on the 
other hand, tq <^ r, i has enough time to interact with 
atom j between two consecutive LREs at site j (atom 
j is frozen on the time scale of atom i). In this case, 
the displacement of atom i affects atom j in the same 
way as in a solid, and the wave propagates. Because the 
vibration of atom i is the superposition of modes with 
different frequencies, we find that modes with frequency 
Lu ^ 1/t {uj <^ 1/t) can (can not) propagate in a liquid. 

We note that this discussion applies to volume- 
conserving shear waves. Low-frequency longitudinal (hy- 
drodynamic) waves propagate in the presence of LREs 
as well, which is related to the existence of the finite 
zero-frequency bulk modulus in liquids (zero-frequency 
shear modulus in liquids is zero by definition). We also 
note that both high-frequency longitudinal and trans- 
verse waves are subject to some damping, but the damp- 
ing of shear waves with frequency uj < l/r is so large 
that it is impossible to talk about their propagation [^. 
Damping decreases with ujt in the ujt > 1 regime for 
longitudinal and shear waves, and increases with ujt in 
the UJT < 1 regime for longitudinal waves [2| . 

In this picture, as we approach a liquid from the solid 
phase, LREs modify the solid vibrational spectrum by ef- 
fectively removing the restoring force (or force constants) 
for shear modes at frequency uj < 1/t. As a result, the 
oscillatory motion starts to "slip" at low frequency, giving 
rise to diffusion due to LREs. On temperature increase, 
the emergent diffusive motion due to LREs contributes to 
the increased kinetic energy of the system K = 3NT/2. 
On the other hand, we can assume that the contribution 
of the emergent diffusive motion to the potential energy 
is small. Lets increase the system temperature T by a 
small amount AT, corresponding to a certain new value 
of T. As restoring forces for shear modes with frequency 
UJ < 1/t are removed, atoms start slipping and diffus- 
ing along the same directions at which they oscillated in 
the presence of restoring forces at temperature T (albeit 
with larger displacements). Note that interatomic inter- 
actions that give rise to the restoring forces at tempera- 
ture T and interactions experienced by the newly diffus- 
ing atoms at temperature T -\- AT have the same micro- 
scopic origin. Therefore, the absence of restoring forces 
at frequency uj < 1/t means that the emergent diffu- 
sive motion weakly contributes to the system's potential 
energy: a large value of this contribution would imply 
strong interactions at frequency uj < l/r and, therefore, 
the existence of restoring forces and shear modes with 
this frequency. Hence, the contribution of the poten- 
tial energy of the emergent diffusive motion to the total 
energy is small compared to other energy terms (the po- 



tential energy of the unmodified longitudinal mode and 
the kinetic energy K), and can be ignored. 



LIQUID ENERGY 



The glass transition temperature Tg is defined from the 



condition r(Tg ) = 10^ s 



Above Tg, each atom par- 



ticipates in the vibrational motion during time r and in 
the diffusional motion when it jumps between two equi- 
librium positions. The instantaneous liquid energy, Ems, 
is the sum of the vibrational, Eyif,, and the diffusional, 
Edif, components: 



Ei, 



E„ 



E, 



dif 



(1) 



Within thermal fluctuations, Eins is equal to the liquid 
energy averaged over time t ^ t. 

Lets now approach the liquid from the solid phase as 
discussed in the previous section. In the solid glass phase 
at T < Tg, Edif = 0, and the energy is purely vibra- 
tional and consists of kinetic and potential energy of the 
oscillatory motion. Lets assume that the temperature is 
instantaneously raised, corresponding to a certain value 
of T < T{Tg) in the equilibrium state. As discussed in 
the previous section, the effect of local rearrangements 
processes is to effectively remove the potential energy 
component that give rise to the restoring force for two 
transverse modes with frequency lu < 1/t. The energy 
of such instantaneously prepared state, E, is the sum of 
the kinetic energy, K, the potential energy of one longi- 
tudinal mode, Ui, and the potential energy of two shear 
modes of frequency uj > 1/t, Ut{i-u > 1/t): 



E^K + Ui + Ut{uJ> 1/t) 



(2) 



We now let the system evolve until the equilibrium 
state is established. In the equilibrium state, E becomes 
the energy of vibrating and diffusing atoms. From the 
energy conservation, the energies of the initial and the 



final state are the same: E = Ei, 



Hence, we can use 



Eq. (2) to calculate the liquid energy. 

Because K = Ki + Kt, where Ki and Kt are the 
kinetic energies of longitudinal and transverse modes, 
E = El + Kt + UtioJ > 1/t), where Ei is the energy 
of the longitudinal mode. If Et is the energy of the 
transverse mode, Kt = Ut — ^ from the equipartition 
theorem. Because the equipartition theorem applies to 
the energy of a single oscillator as well as to the sum 
over any set of oscillators, C/t(w > 1/t) — *'■"> ''^i ^ 



Then, E = El + ^ + E,(u>i/r) 



Et can be writ- 



2 ' 2 • -^t 

ten as the sum over modes below and above uj — 1/t, 
Et = Et{uj < 1/t) + Et{LJ > 1/t), where the two terms 
refer to their solid-state values. Then. 



, , Etlcj < 1/t) , , 

E = Ei+ Et{uj > 1/t) + ^ (3) 

The first two terms in Eq. (3) are the vibrational ener- 
gies of one longitudinal mode and two transverse modes 
with frequency uj > 1/t. In the solid glass phase, where 
T is very large, these two terms become the vibrational 
energy of the system, and the last term is zero. On in- 
creasing the temperature in a liquid, the second term 
decreases. The last term ensures that the total kinetic 
energy of the system does not change as a result of this 
decrease: the kinetic energy is defined by temperature 
only {K = 3iVT/2) regardless of how it partitions into 
vibrational and diffusional motion. 

To calculate the first two terms in Eq. (3), we sep- 
arate N longitudinal and 2N transverse normal modes 
with frequency uj > 1/t in the statistical sum of the vi- 
brational motion: 



2 2\ 

- '^uQi ) 



dpdq (4) 



f 1 ^ 

X / exp I -— J2 (^? + ^''«^') ) '^P'^1 

where ujn and loh are frequencies of longitudinal and 
transverse vibrations and loq = 1/t. This gives 



z = Zi ■ ZtiLu > 1/t) = (2^r)^ n ^i 



N 



, where Zi and Zt{uj > 1/t) are 



2N 

(2^r)^i ( n ^u 

the contributions to Z due to the longitudinal mode 
and the transverse modes with frequency lo > 1/t, re- 
spectively, and A^i is the number of transverse modes 
with LO > loq. The corresponding energies are Ei — 
T^-^ln{Zi) ^ NT and Et{Lu > 1/t) = T'^-^lniZt{uj > 
1/t)) = NiT. This gives the first two terms in Eq. (3). 
iVi can be calculated from the density of states in 
the Debye model. Recall that in an isotropic solid 
(glass), vibrational density of states is written as g{u!) — 

uP' , where c/ and Ct are longitudinal and 
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transverse sound velocities [l| . The Debye frequency ujm 
is defined from the condition J g{uj)duj — 3N, where 



N is the number of atoms in the system. If the mean 

\, the density 



2 



speed of sound is introduced as -^ 

of states becomes g{Lu) = ^^" , where a;f„ — 6tt^Nc'^/V 
1]. Instead of a;„i, lets define transverse Debye frequency, 

LUmt, using the condition J gt{uj)duj — 2N, which re- 



fiects the fact that 2A^ degrees of freedom are associated 
with transverse vibrations. Here, gt{uj) — -^-^^"^ is 



the density of states of transverse vibrations. This gives 
^3 = Qtt'^Nc^JV and gAoj) = -^w^. Then, A^i = 

J' gt{Lj)dcj = 2N (l- (^) J, and the vibrational 

energy of one longitudinal and two transverse modes with 
frequency lu > 1/t becomes Eyu, ~ Ei + Et{uj > 1/t) = 

{N + Ni)T = TVT ( 3 - 2 i-^X 

To calculate the last term in Eq. (3), we note that 
similarly to Et{uj > 1/t) = NiT, Et{uj < 1/t) = N2T, 
where N2 is the number of transverse modes with uj < 

iVQ. Because N2 ^ 2N - Ni, N2 = 2N (^] . Then, 
E:^{N + Ni+ N2/2)T, giving 



HEAT CAPACITY 



E ^ NT \3 




\ ^mt 



(5) 



Eq. ([5]) predicts the relationship between the energy 
of a liquid and its relaxation time r — l/t^o- The first 
term in Eq. ([5]) is the classical energy of a solid. The 
second term is the contribution to the energy due to lo- 
cal rearrangement processes in a liquid. According to 
Eq. ([5|), liquid energy is equal to the solid-state value of 
iNT in the very broad range of relaxation time t. This 
range starts from the experimentally accessible r = 10'^ 
s, and continues up to the Debye vibrational period of 
td = 10^^'^ s. Only when cjq becomes of the order of 
Debye frequency, liquid energy and heat capacity start 
decreasing from their solid-state values. Hence, Eq. ([5]) 
predicts that this decrease should start in low- viscous liq- 
uids not far above their melting point. This agrees with 
experiments on liquid Hg, K, Na and other systems [8|,|9|. 
If L0{) = (jJmt ~ 1/''D, Eq. ([5]) gives the lower limit of the 
energy of a simple monatomic liquid of 2NT . This is the 
result that Brillouin obtained under t he assumption of 
the complete loss of transverse modes 
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The reason why the energy of a liquid is equal to the 
solid-state value of iNT unless o^o is comparable to the 
Debye frequency can also be discussed on the basis of 
Eq. (1). At any given moment of time, the number of 
particles in the transient diffusional state, Ndif, is N^if = 
A'o exp(— y/T), where V is an activation barrier for a 
LRE {V can be temperature-dependent) and Nq is the 
total number of atoms. Because r = td exp{V/T) (see 
Ref. i), ^ = m. When r » tq, Nd^f/No « 1, 
i.e. most of the particles in the system are vibrating, 
and only a negligibly small proportion of particles are 
diffusing. Because Ndif/No <C 1, the contribution to the 
energy from diffusion can be ignored in Eq. (1), giving 
Eins — Evib- As we have seen above, the vibrational 
energy, E^it = {N + Ni)T = 3NT when t > rt), i.e. the 
same as in the solid. This point will be discussed below 
for its implications for the problem of glass transition. 



From Eq. ©, c„ = ;^^ is: 



dT 



T 3 



Tmt(jo 




(6) 



where Tmt is the Debye period for transverse vibrations, 
T,nt ~ TD, 77 is viscosity, and we used the Maxwell re- 
lation 77 — tGoo {Goo is instantaneous shear modulus) 
in order to directly relate c„ and viscosity. Note that 
Eq. ^ relates heat capacity and viscosity with no fit- 
ting parameters, because T^tG 00 is given by the liquid 
properties. 

Mercury is an appropriate liquid to test the predic- 
tion of Eq. ^ for a number of reasons: it is a sim- 
ple monatomic liquid, anharmonic effects in mercury are 
small [9| (note that thermal expansion is ignored in our 
discussion), and finally, it is a low- viscous liquid, which 
makes it possible to observe the decrease of heat capacity 
in the reasonable temperature range: as discussed above, 
this is expected to take place when loq becomes of the 
order of Debye frequency. 

We have taken measurements of c^ of mercury from 
Ref. [8|, which has electronic contribution subtracted, 
hence giving heat capacity due to ions only. Viscosity 
data was taken from Ref. [20|, and interpolated in or- 
der to use in Eq. © to calculate c^. Because viscosity 
shows the deviation from the Arrhenius behaviour, we 
used the Vogel-Fulcher-Tammann expression to fit the 
data, 77 — 77oexp(A/(T — To))- A good fit was achieved 
with 770=0.612-10-3 Pa-s, A=235.157 K and To=40.726 
K (see Fig. la). Because the experimental temperature 
range of viscosity is 250-600 K [20|, slightly smaller than 
the experimental temperature range of c^ , 234-754 K [8[ , 
we used the fitted expression to extrapolate viscosity data 
in the temperature range of c^ (see Fig. 1). 

Eq. ([6]) has no fitting parameters, hence there is no 
fiexibility to fit its prediction to the experimental curve 
of Cy. However, because Tmt and Goo are not known pre- 
cisely for mercury, in practice we can use T^tG 00 as a 
fitting parameter. In Fig. lb, we compare c„, calcu- 
lated from Eq. ([6]), with the experimental data, using 
T„tGoo=0.55-10-3 and 0.57-10-3 Pa-s. Given that the 
experimental error of Cv is 0.1-0.2 J/K [8| (experimental 
error of viscosity is not known), we find a good agreement 
between the calculated and experimental Cy . 

It is important to note that the used value of TmiGoo in 
Fig. lb gives physically sensible values for Tmt and Goo: 
if Goo of mercury is about 5 GPa [21], the used values of 
TmtGoo imply Tmt of about 0. 1 ps, a typical value of Debye 
vibrational period. We also note that the used value of 
TmtGoo is close to that in other liquids: in sodium, for 
example, the product of high-temperature r and Goo is 
about 0.5-10-3 Pa-s 141. We can therefore conclude that 
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FIG. 1: (a) Circles are experimental data of viscosity of mer- 
cury, line is the fit and extrapolation using the Vogel-Fulcher- 
Tammann formula; (b) Solid line is the experimental c„ for 
mercury. Dashed lines are calculated values, using TmtG oo of 
0.55 -10"^ Pa-s (top) and 0.57 -10"^ Pa-s (bottom). 



in the previous sections is at odds with this assump- 
tion. Indeed, we have discussed that transverse waves 
with frequency a; > l/r propagate in liquids, and their 
frequency range only increases at low temperature. We 
have seen that transverse waves disappear at high tem- 
perature only, whereas neglecting them at low tempera- 
ture gives, in classical case, C^ = "^N, inconsistent with 
the experiments. Interestingly, we show that transverse 
waves can not be excited in a quantum liquid in the low- 
temperature limit, confirming the Landau assumption. 

In a quantum solid, Eyn, is the sum of energies of quan- 
tum oscillators, represented by one longitudinal and two 
transverse modes. In a quantum liquid, Eyn, is the sum 
of energies of one longitudinal and two transverse modes 
with frequency lo > l/r. As in the Debye model of a 
solid, the summation over the phonon frequencies can 
be substituted by integration using the Debye density 
of states [1[ . Similar to uJmt introduced previously in the 
density of states of transverse vibrations, gt{^) = -t^o;^, 
we define the Debye longitudinal frequency, w„i/: from 



the condition 



1 



N, which reflects the 



fact that N degrees of freedom are associated with lon- 
gitudinal vibrations. This gives w^; = 6^T^Nc^/V, and 
the density of states of longitudinal vibrations, gi^to) — 
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Evib becomes 
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(7) 
where Eq is the energy of zero-point vibrations and 
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the decrease of Cy follows, within the experimental error, 
the decrease of viscosity reasonably well, with physically 
realistic values of Tmt and Goo ■ 
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VIBRATIONAL ENERGY OF A QUANTUM 
LIQUID 

In most liquids, crystallization or vitrification take 
place at low temperature. However, if a liquid state exists 
at low temperature, quantum effects become important. 
Below we calculate the effect of relaxation process on the 
vibrational energy of a quantum liquid, Eyij,. 

We first note that in considering the phonon exci- 
tations in a quantum liquid at very low temperature. 
Landau assumed that these are due to the longitudinal 
phonons only [Ij. The vibrational energy is then three 
times smaller than that of a quantum solid, reflecting the 
assumption that excitations due to two transverse modes 
are completely lost [II]. At first glance, the discussion 



where D{x) 



■I 



A^ 



— is the Debye function [l|. 



The first three terms in Eq. (8) is the vibrational en- 
ergy of a quantum solid, obtained from Eq. ^ ii ujq — 0. 
The last term in Eq. (8) is the contribution from the 
relaxation process in a liquid. In the high-tcmpcrature 
limit, when D{x) — 1, E^n, becomes the vibrational en- 
ergy of a classical liquid that follows from Eq. (4). 

The low-temperature limit of the Debye model corre- 
sponds to the large value of the argument oi D{x) , so that 

4 

neglecting exponentially small terms, D{x) — -^ and 
E„f, ex T" [l|. In Eq. (8), this takes place when ^ > 1 

(which also imphes ^"^' :s> 1, because luq < uJmt)- When 
D{x) oc T-, the last two terms in Eq. (8) disappear for 
any loq, corresponding to the complete loss of excitations 



due to transverse waves. This can be understood because 
at temperature T, only phonons with frequency fiiij < T 
are excited and contribute to E^ and the temperature- 
dependent energy terms in Eq. ([7]). Hence, in the low- 
temperature limit Hloq ^ T only phonons with frequency 
uj <S^ loq are excited, but transverse phonons with this 
frequency are excluded from Eq. ([7]). As a result, trans- 
verse modes can not be excited in a liquid in the low- 
temperature limit, and all excitations in the system are 
due to the longitudinal phonons only. 



IMPLICATIONS FOR GLASS TRANSITION 

There are two implications of the discussed picture for 
the long-standing problem of glass transition 19|, |23, [24 1 . 
The first one concerns the general approach to glass tran- 
sition from a liquid state. Basic discussions of liquids be- 
gin with stochastic dynamics, which describes Brownian 
motion and self-diffusion, and is modeled by Langevin- 
type and related equations [J, la, LZ|- These approaches 
assume stochastic dynamics, corresponding to the regime 
r ~ td. To discuss glass transition, the formalism is then 
modified so that is can be extended to the supercooled 
state where r 3> td and eventually to the glass phase 
where riTg) = 10^ s. The important point here is that 
the condition r ^ td holds in 15 decades of r (10~^^-10^ 
s). In this range of t, the system is a liquid as follows 
from its inability to support static shear stress on exper- 
imental time scale, i.e. it flows, yet its main properties 
are those of a solid. First, as discussed above, the hquid 
energy is equal to the vibrational energy of a solid in the 
range of 15 decades of r. Second, a hquid supports sohd- 
like shear modes with frequencies oj > 1/t which can 
also span up to 15 orders of magnitude depending on r. 
Hence, in almost entire glass transformation range, the 
important properties of a liquid are governed by its solid- 
like elastic features. This suggests that a glass transition 
theory can include the solid-like elastic properties of a 
liquid from the outset, rather than modify and extend 
the stochastic high-temperature behaviour (where these 
properties are absent) to the supercooled liquid state and 
glass. Basing on this view, we have recently proposed 
how glass transition can be understood on the basis of 
liquid elasticity [25|. 

The second implication concerns the behaviour of liq- 
uid heat capacity at Tg. Opposite to the increase of c„ 
on lowering the temperature in Fig. 1, constant-pressure 
heat capacity decreases at Tg, often with a jump, Acp 
19l l23l . |24J . It is therefore important to discuss how the 
two effects can be reconciled. Several popular theories of 
glass transition relate Acp at Tg to the reduction of liq- 
uid configurational entropy and to an underlying phase 
transition [l9| . This approach has been convincingly crit- 
icized for a number of reasons [19| . We propose that Acp 
at Tg can be understood as a kinetic effect related to the 



non-equilibrium nature of glass transition, rather than to 
the presence of an underlying phase transition. Hence, 
this effect is different from the one that governs liquid 
heat capacity in Eq. (6): recall that we related liquid 
heat capacity to the liquid equilibrium vibrational states. 

Acp at Tg can be understood from the disappearance 
of LREs at Tg. The presence of these processes in a liq- 
uid (or their absence in a glass) has a profound effect 
on several of important properties of the system. First, 
LREs are accompanied by the increases of local volume. 
This is because the energy needed for an atom to es- 
cape its cage at the constant volume is very large due to 
strong short-range interatomic repulsions. Hence, atoms 
in the cage need to increase its volume in order to allow 
for the escape of the central atom [2, |l9| . Because Tg is 
defined from the condition of disappearance of LREs on 
the experimental time scale, thermal expansion is larger 
(smaller) above (below) Tg. This gives rise to the exper- 
imentally observed kink of density at Tg [26|, and con- 
tributes to the jump of Cp at Tg. Second, because LREs 
are accompanied by fluctuating local volume increases, 
they result in a faster softening of the high-frequency 
elastic constants of a liquid as compared to a glass, as 
witnessed by the kink of the high-freq uency longitudinal 
and transverse sound velocities [271, l28| . This introduces 
a difference in the compressibility between a liquid and 
a glass at Tg, and contributes to Acp. Finally, at times 
t ^ T, LREs ensure that the configurational entropy of 
the system, Sc, is equal to the "entropy of melting" asso- 
ciated with the increase of the number of configurations 
of a liquid over that of a solid [3] . If LREs are frozen at 
Tg, the change of Sc contributes to Acp. 

Importantly, small variations of temperature around 
Tg give large changes of r, leading to the freezing of 
LREs in a very narrow temperature range and therefore 
giving rise to the observed apparent jump of Cp. Gen- 
erally, this is because at Tg, the activation barrier for 
LREs, V, is very large: because r = td exp{V/T) and 
T{Tg) =100-1000 s, V/Tg =35-37. Moreover, in "frag- 
ile" liquids [23[, i.e. liquids with large slope of ln(r) vs 
Tg/T at Tg, small variations of temperature around Tg, 
AT, result in a large change of t. Take, for example, 
propylene carbonate, for which relaxation time follows 
the VFT law, r = td exp(^/(T — Tq)) with parameters 
A = 612 K, To = 139 K and log(TD) = -12.3 s ^. The 
conditions T{Tg) = 100 s and T{Tg - AT) = 1000 s give 
Tg « 157.6 K and AT « 1.2 K, i.e. the decrease of tem- 
perature at Tg by only 1.2 K gives a ten- fold increase of 
T. Similarly, for several other fragile liquids (29|, we find 
that the variations of AT — ±(1-2) K at Tg result in the 
variations of t by an order of magnitude. 

This picture predicts that no Acp would have been ob- 
served had the experiment lasted for time t ^ t. This 
is consistent with the general finding that the observed 
anomalies at Tg depend on cooling rate and observation 
time. Another prediction of this picture is that Acp at 



Tg should be sharp in fragile, but not in "strong" liquids, 
which have a smaller slope of ln(T) vs Tg/T at Tg. This 
is consistent with the experimental results [23| . To sum- 
marize, Acp can be understood as a natural signature of 
the glass transition insofar as this transition is defined by 
the freezing of LREs at the experimental time scale, but 
not as a result of a phase transition. 



CONCLUSIONS 

In summary, we discussed the solid-like elastic prop- 
erties of liquids, including its vibrational states, energy 
and heat capacity. The decrease of liquid heat capacity 
with temperature was attributed to the increasing loss of 
two transverse modes with frequency lu < 1/r. In a sim- 
ple model, we related liquid heat capacity and viscosity, 
and compared this relation with the experimental data 
of mercury. We also calculated the vibrational energy of 
a quantum liquid, and showed that transverse phonons 
can not be excited in the low-temperature limit. Finally, 
we discussed the implications of the proposed approach 
to liquids for the problem of glass transition. 
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